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Abstract: In view of the recent excitement about a tension between determinations of 
i^Ds from experiment and from simulations of lattice QCD with dynamical quarks, we try 
to clear up the picture of lattice determinations in the continuum limit of the quenched ap- 
proximation. For 0(a) improved Wilson quarks we see linear scaling in the squared lattice 
spacing a? only for a < 0.08 fm. For coarser lattices we observe significant contaminations 
from higher order cutoff effects. As an aside we also study the scaling of the charm quark 
mass and the ratio of the vector to the pseudo-scalar decay constant and the spin-splitting. 
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1. Introduction 

Since the time when simulations of lattice QCD with dynamical fermions became feasible, 
the phenomenology of the Dg-meson has been considered as a field where lattice QCD 
could provide benchmark predictions that would be confronted with increasingly precise 
measurements from experiments like Belle, BaBar and CLEO-c. Preparatory tests of the 
techniques in quenched lattice QCD [1,2] indicated that a precision of only a few percent 
for observables like Fd^ is possible when simulating the full theory, while keeping all sources 
of systematic uncertainties under control^ . 

Only recently, precise determinations of F^^ appeared in dynamical simulations with 
Nf = 2 and 2+1 flavours of sea quarks [3-6]. The central values, although partly pre- 
liminary [5,6], turned out to be about 10% larger than in the quenched case. This was 
not surprising given that the quenching error was always estimated to lie in the range of 
10-20%. Surprising, however, turned out to be a recent comparison with a compilation 
of experimental results for F^^ [7] by CLEO-c [8] and Bcllc [9, 10]. A tension between 
the experimental value and lattice results triggered speculations about possible signs for 
physics beyond the Standard Model [11]. 

In this work we want to emphasize that the lattice spacing dependence must be mapped 
out over a large range in order to unambiguously isolate the leading lattice artefacts — 
this is all the more important in the heavy quark sector. In particular the approach 
to the continuum limit of the leptonic decay constant of the Dg-meson in the quenched 

^The systematic effect due to the chiral extrapolation of the sea quark masses cannot be assessed in the 
quenched theory. 
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approximation, as recently summarized in [2,12], did not exhibit proper scaling, and we 
try to clarify this issue here. It appears that simulations with lattice spacings larger than 
a 0.08 fm can suffer from ambiguities in the 0(a) improvement of the quark bilinear 
currents, which significantly affects the extrapolation of the data to the continuum limit. 
It is indispensable to go to smaller lattice spacings in order to assure scaling. These findings 
are extremely important for assessing current error estimates in lattice studies of Fd^ in the 
unquenched case as they have been reported by ALPHA, ETMC, HPQCD and MILC [3-6]. 

To this end wc extend the ALPHA Collaboration's computation of F^^ in the quenched 
approximation reported in [1] for lattice spacings in the range a 0.09 — 0.05 fm by sim- 
ulations of an additional lattice spacing a Ri 0.03 fm [13]. Moreover, we increased the 
statistics with respect to [1] and [13] significantly, which allows us to present a comprehen- 
sive picture of the lattice spacing dependence of the decay constant i^Os in quenched lattice 
QCD with non-perturbatively improved Wilson quarks. In addition, we also present scaling 
studies for the renormalization group invariant charm quark mass Mc [13, 14], the mass 
splitting mDj — and of the ratio of the vector to pseudo-scalar meson decay constants 
Fd*/Fd^ [13]. We start by briefly explaining the setup of our large- volume simulations of 
the quenched QCD Schrodinger functional [15, 16] and then discuss the simulation param- 
eters. This is succeeded by a presentation of our results, before ending with a discussion 
and our conclusion. 

2. Details of the simulation 

In the following we introduce the set of Schrodinger functional correlation functions, from 
which we determine niD^ , mj)* , Fd^ , i<b* and the mass of the charm quark. We compute the 
correlation functions directly for all required mass-degenerate and non-degenerate combi- 
nations built from quark propagators at the physical strange quark mass and at two values 
of the charm quark mass close to the physical one. This allows for an interpolation of the 
results to the exact physical point, which we define through the experimentally determined 
mass of the Dg-meson as input. This computation is repeated for five different lattice 
spacings in the range of a ~ 0.1 — 0.0 3 fm in order to investigate the scaling behaviour of 
our observables and identify a range in the lattice spacing where observables scale linearly 
in o^. 

2.1 Correlation functions 

For this study the relevant Schrodinger functional heavy-light 2-point correlation functions 
are 

fA{xo) = -UM^)0) , /p(xo) = -^(P(x)0), /i = -^(C»'0), 

fy{xo) = -lUVi{x)Oi) , hA^o) = -lUTi^iOiix)) , k, = -^j:{Opi), ^^-^^ 

i i i 

where Af^(x) = '4' s{x)l iilb'4'hix) is the axial-vector current, P{x) = V's(^)75V'h(^) is the 
pseudo-scalar density, and V^(x) = V's(^)7M^h(^) a-i^d Tfj_y{x) = ^six)'^fiu'^h{x) are the 
vector and tensor currents. In each case we indicate the respective quark flavour with a 
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subscript "s" for light (= strange) and "h" for heavy quarks with masses close to the one 
of the charm quark. The quark bilinears 

^ = 1^ E Ch(y)75Cs(z) , o' = ^J2t{yHCL{^) , (2.2) 
= ^ECh(y)7iCs(z), oi = |jEcl(y)7^Ch(z), (2.3) 

are the meson wall sources at the xo = and the xq = T boundary timeslices, respectively. 
For more details on Schrodinger functional correlation functions and for unexplained no- 
tation we refer the reader to, e.g., refs. [17,18]. We employ the improved axial- vector and 
vector currents 

Al{x) = Aq{x) + CAa5oP(x) and F/(x) = Vi{x) + cvadoToi{x) , (2.4) 

where Bq is the time component of the symmetrized next-neighbour lattice derivative. The 
improvement coefficients ca(5o) cind cy{go) for the quenched theory have been determined 
non-perturbatively in [19] and [20,21], respectively^. On the lattice, the axial- vector current 
and the vector current receive a multiplicative (scale independent) renormalization and we 
write 

{An)o = Za[1 + 6A(amq,h + amq,s)/2] 4 + 0(0^) , (2.5) 
{VR)i = Zy[l + bviam^^i, + amq,s)/2] + O(a') . (2.6) 

The corresponding renormalization constants for the quenched theory ZA{go) and Zy{go) 
have been determined non-perturbatively in [23]. 0(a) artifacts that are proportional to the 
bare subtracted quark masses amq j = ^(^ — ^r^) for i = s,h, where Kent is the critical 
hopping parameter, are canceled by terms proportional to the improvement coefficients 
bA{go) and bv(5'o); they have been calculated in 1-loop perturbation theory [20]. 

2.2 Observables 

We define the pseudo-scalar effective mass via 

ampsixo + a/2) = log ' (^.7) 

fX{xo + a) 

and the expression for the local pseudo-scalar meson decay constant is given by [18] 

aFpsixo) = -2Za [1 + 6A(amq,h + amq,s)/2] (mpgL^)"^/' e(-°-'^/2)-Ps 

v/i 

X |l - r/X^e-^o^ - j?^e-(^-^°)"*G} + 0(0^) , (2.8) 

and an analogous expression holds for the mass and the decay constant in the vector 
channel. In (2.8) the factor (mpsL^)~^^'^ takes into account the normalization of one- 
particle states, and f^^^"^ cancels out the dependence on the meson sources. Because of 



^In the case of cv(flo), the parameterization derived in [22] has been used here. 
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its exponential decay, the correlation function fj^ is dominated by the ground state in the 
pseudo-scalar channel for large time separations from the boundaries (i.e. xq and T — xq). 
Hence, eqs. (2.7) and (2.8) are expected to exhibit a plateau at intermediate times, when 
the contribution ri^e~^°^ of the first excited state and the contribution ri'j^Q-{T-xo)mG 
from the 0++ glueball both are small and well below the statistical errors, such that final 
estimates for masses and decay constants are obtained by fits to a constant within the 
plateau regions of the local quantities (2.7) and (2.8), as explained in [18]. 

Concerning the quark masses, we employ the continuum PCAC relation, which in the 
situation of QCD with non-degenerate quarks at hand can be written as 

d^A^{x) = 2mhsP(x) = (mhh + mss)P{x) (2.9) 

and on the lattice translates into a definition of the bare current quark mass in terms of 
the heavy-light correlation functions (2.1) [17], 



arribsixo) = | 



^{do + do)fA{xo) + CAadodofp{xo) 



//p(xo) , (2.10) 



do (9q) being the forward (backward) lattice derivative in the time direction. In this 
notation, a single quark mass ma, i = s,h, is obtained if the quark flavours are mass 
degenerate. Again it is understood that the bare heavy-light current quark masses entering 
the formulae below are computed as timeslice averages within a central flat region of the 
associated local masses, eq. (2.10). 

After multiplicative renormalization of the improved axial current and the pseudo- 
scalar density we obtain the identity 

r , N „ -^A [l + &A|(amq,h + amq,s)] . r^, 2^ 

a(mR,h + mR,s) = 2 . , ■ -f arrihs + 0(a ) 

Zp [1 + 6p2(amq,h + amq,s)J ^2 n) 

= 2ZaZp^ [l + (6a - bp)l(amc^,h + aniq^s)] anihs + O(a^) 

for the sum of the renormalized PCAC heavy and strange (valence) quark masses in the 
0(a) improved theory, where (b^ — bp){go) has been determined non-perturbatively in [24]. 
The scale dependent (and, in quark mass independent renormalization schemes, also flavour 
independent) renormalization factor Zp = Zp{go, L/a) encodes the scale dependence of the 
renormalized quark masses and is available from [25] in the Schrodinger functional scheme 
for A^f = and the range of bare couplings relevant here. Eventually, using the flavour 
independent ratio M/m-n^idi), /i = 1/L, in the continuum limit, we combine this with 
Zp^jZp into the total renormalization factor .^m(5o) [25] to express our results directly in 
terms of the scale and scheme independent renormalization group invariant (RGI) quark 
mass M. The RGI mass is very convenient since it may then be straightforwardly converted 
to any other renormalization scheme, such as MS at some desired scale, by means of 
continuum perturbation theory. 

From the various quark mass definitions for non-degenerate and degenerate quarks we 
construct the following expressions for the 0(a) improved RGI heavy (i.e. charm) quark 
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mass. Once we compute it via the current quark mass from the heavy-strange current with 
the strange current quark mass subtracted, 



- aniss [1 + (6a - 6p)amq,s]| , 



(2.12) 



and once from the case of degenerate quarks via the current quark mass from the heavy- 
heavy current, 

aMj^^^ = Zm amhh [1 + (6a - 6p)amq_h] • (2-13) 

A third definition of the RGI charm quark mass is given directly in terms of the bare 
subtracted heavy quark mass through the relation 

aMj^^ = ZmZ amq,h [1 + 6inamq,h] • (2.14) 

We will use the quenched parameterization of ZM{go) from [13], which extends the one 
based on [25,26] to also include (3 = 6.7859 (see below), and of Z{gQ) = {Zy^Z-p / Z jCjigo) 
and 6111(50) from [24]. These three definitions of the quark mass will have different cutoff 
effects, but should agree with each other once being extrapolated to the continuum limit. 
This provides a nice check of the extrapolation procedure. 

2.3 Parameters for the scaling study 

We have generated two ensembles of gauge field configurations (called A and B in the 
following) with standard algorithms, for five different lattices of approximately constant 
physical size L/vq ~ 3 in spatial directions, but decreasing lattice spacing. Here, tq is the 
Sommer scale [27] which we use in order to estimate the lattice spacing a in physical units. 
Some lattice and simulation parameters arc collected in table 1. The numerical simulations'^ 
were carried out on the APEmille and apeNEXT computers at DESY Zeuthen and on the 
IBM p690 computers of the HLRN [28]. In the latter case we have employed an adapted 
and performance-improved version of the Schrodinger functional implementation based on 
the MILC code [13,29]. 

2.4 Hopping parameters 

The values of the critical hopping parameter Kcrit, which enter the analysis through the 
quark mass dependent 0(o) terms in eqs. (2.5) and (2.6) and the subsequent expressions, 
were gathered from the non-perturbative determination in quenched QCD of [20] and, 
where necessary, we interpolated Kcrit to the desired value of the lattice coupling [13]. 
Kg was fixed prior to the simulations using published results for the PCAC relation 

M, + M = ZM^mi (2.15) 



*A11 simulations were done in double precision arithmetic, except for the subset of simulation points of 

data set A which wore the basis of [1]. However, we did not notice any visible impact of the precision used 
for the charm quark propagator computation on our results. 
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data set A 






data set B 








L/a 


T/a 


L/ro 


a [fm] 




^update 


noR 




'^update 




Pi 


= 6.0 


16 


32 


2.98 


0.093 


380 


100 


8 


2100 


25 


5 




= 6.1 


24 


40 


3.79 


0.079 


201 


100 


12 


1300 


25 


5 




= 6.2 


24 


48 


3.26 


0.068 


251 


100 


12 


1300 


25 


5 


Pa 


= 6.45 


32 


64 


3.06 


0.048 


289 


100 


12 


1400 


40 


10 


P5 


= 6.7859 


48 


96 


3.00 


0.031 


150 


50 


24 


604 


50 


20 



Table 1: Lattice geometries and simulation parameters used in the simulations for the scaling 
study for data set A and data set B. The correlation functions were evaluated on rimeas gauge 
configurations, which were separated by riupdate update iterations consisting of one heat-bath and 
noR over-relaxation sweeps. 



for 0(a) improved Wilson fermions in quenched QCD [30], where Mg is the RGI quark 
mass of the strange quark and M = |(Mu + M^) the average RGI light quark mass. Fk is 
the kaon decay constant, and Gk denotes the vacuum-to-K matrix element of the pseudo- 
scalar density. Note that in contrast to [1] we here have determined Kg for each value of 
(3 by fixing the RGI strange quark mass to its value found at finite lattice spacing in [30] . 
Therefore, differs from the one in [1] by 0{a^). Wc take over the values for Zu ^ 
from [30] at /3 G {6,6.1,6.2,6.45} and extrapolate those linearly, in order to arrive at the 
estimate 1/3=6.7859 = 0.2268(57). Using (2.15) and the ratio Mg/M = 24.4±1.5 from 

continuum chiral perturbation theory [31], the values for Ks{P) are then determined as the 
solutions of 

ro(Ms + M) = ZmZ amq,s [1 + 6mamq,s] (l + , (2-16) 

where the r.h.s. is meant to be read as a function of /3 = G/qq- 

The choice of hopping parameters of the charm quark, , relics on the results for Mc 
at P G {6,6.1,6.2,6.45} in [14], which in addition were extrapolated linearly in (a/ro)^ to 
P = 6.7859. The estimates for in table 2 have then been obtained by solving a quadratic 
equation similar to (2.16). We have also generated data for the supplementary hopping 
parameters and (for data set A and data set B, respectively), which yield results 
close to • This enabled us to interpolate all observables to the point corresponding to 
the physical charm quark. 

3. Analysis and results 

We analyzed our data using the T-method [32] , where the statistical errors of the observ- 
ables are estimated by directly analyzing autocorrelation functions, and cross-checked its 

outcome by a jackknife procedure. Generically, autocorrelation times turned out to be 
small so that our measurements on gauge field ensembles at given (5 could be treated as 
statistically independent. 
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(3 


6 


6.1 


6.2 


6.45 


6.7859 


/^crit 


0.135196 


0.135496 


0.135795 


0.135701 


0.135120 


Kg 


0.134108 


0.134548 


0.134959 


0.135124 


0.134739 




0.119053 
0.128790 


0.122490 
0.130750 


0.124637 
0.131510 


0.128131 
0.132690 


0.130823 
0.132440 


I^C2 


0.123010 


0.125870 


0.127470 


0.130030 


0.130253 



Table 2: Summary of all hopping parameters. 



/5 


fi) 

'•()»' i>s 


data set A 

(2) (1) 


(2) 

'•o'».v 


data set B 

(1) (2) 


6.0 


4.295(28) 


4.970(33) 


4.644(32) 


5.255(36) 


4.293(27) 


4.971(31) 


6.1 


4.273(29) 


4.987(34) 


4.618(33) 


5.271(37) 


4.277(29) 


4.993(33) 


6.2 


4.286(31) 


5.013(37) 


4.651(37) 


5.316(41) 


4.282(31) 


5.012(36) 


6.45 


4.282(36) 


5.026(42) 


4.655(42) 


5.332(46) 


4.261(35) 


5.001(41) 


6.7859 


3.738(37) 


5.182(51) 


4.139(46) 


5.462(55) 


4.825(46) 


5.181(50) 



Table 3: Results for the pseudo-scalar and vector meson masses for the data sets A and B. 
Superscripts "(1)" and "(2)" refer to hopping parameters and (see text). 





A 


A+B 


A 


A+B A+B A+B 


(3 


roirriB* - mo J 






roMi'^^) roMi^^) roikfi^) 


6.0 


0.280(14) 


0.5165(62) 


1.018(42) 


4.369(47) 5.203(56) 3.236(35) 


6.1 


0.285(11) 


0.5604(67) 


1.014(39) 


4.272(46) 4.803(52) 3.477(38) 


6.2 


0.314(15) 


0.5809(84) 


1.131(68) 


4.243(47) 4.640(51) 3.689(41) 


6.45 


0.324(17) 


0.5813(76) 


1.142(49) 


4.170(47) 4.368(50) 3.931(45) 


6.7859 


0.341(13) 


0.5613(81) 


1.064(50) 


4.070(48) 4.171(50) 3.980(48) 


c.l. 


0.346(15) 


0.557(11) 


1.066(61) 


4.040(56) 4.055(58) 4.090(54) 



Table 4: Results for the mass splitting, decay constants and the RGI charm quark mass at finite 
lattice spacing and in the continuum limit (c.l., bottom row). 



The results for the pseudo-scalar and vector meson masses are summarized in table 3. 
For both of the data sets A and B we have computed results for pseudo-scalar masses very 
close to the physical point romog = 0.5 fm x 1969 MeV ^ 4.988 [33]. In order to make 
predictions at the physical point, we interpolate the mass splitting and the RGI charm 
quark mass linearly in the pseudo-scalar meson mass and the decay constant linearly in 
the inverse pseudo-scalar meson mass, as suggested by HQET. In the pseudo-scalar channel 
we have results for primary observables from the (statistically independent) data sets A 
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0.55 
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° this work ca by ALPHA [20] 
this work ca by LANL [21] 
Ah-Khan et al. [2] 




I I I I I I I 







0.03 



0.01 0.02 

Figure 1: Extrapolation to the continuum hmit of the decay constant F^^. We also show the 
results of the recent simulation with non-perturbativcly improved Wilson fermions in [2] . Only the 
data represented by filled squares entered the fit. 
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□ CA by ALPHA [20] 
ca by LANL [21] 



1-2 



1.1^ 
1 

0.9 




I I I I I I I - 



0.01 0.02 0.03 

(a/ro)2 

Figure 2: Extrapolation to the continuum limit of the ratio Fd» /-Fbs- Only the data represented 
by filled squares entered the fit. Note that this quantity was extracted from data set A only, as the 
vector channel correlators are available only in this case. 



qO.35 



o 




0.25 



0.01 0.02 0.03 

Figure 3: Extrapolation to the contimmm limit of the mass splitting. Only the data represented 
by filled squares entered the fit. As in figure 2, this quantity was extracted from data set A only. 



and B; therefore, we average the results and add the statistical errors in quadrature (cf. 
table 4). 

4. Discussion 

The approach of the data to the continuum limit for toFds, -Pb|/-^Ds) '"o("iD* — "^-dJ 
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0.01 0.02 0.03 



(a/ro)2 

Figure 4: Extrapolation to the continuum limit of the RGI mass of the charm quark. Only the 
data represented by filled symbols entered the fit. For the red star and squares, the ca of [20] 
(ALPHA) was used, while it is the one of [21] (LANL) for the grey data points. 

and the three definitions of the RGI charm quark mass is depicted in figures 1-4 for all 
five lattice spacings. As far as they enter, both sets of non-perturbative values for the 
improvement coefficient ca(5'o) of the axial- vector current, i.e. from ref. [20] (ALPHA) and 
from ref. [21] (LANL) and differing in the particular improvement conditions imposed for 
their determination, have been considered. 

In ref. [35] a breakdown of 0(a) improvement for improved Wilson quarks with heavy 
quark masses about aAfh > 0.64 was found in perturbation theory. Since the charm 
quark mass in lattice units is heavier than this bound for the two coarsest lattices with 
a = 0.093 fm and a = 0.079 fm, we exclude the results on these lattices from the continuum 
extrapolations. In fact, significant deviations from the linear dependence on a? of the data 
on the coarser lattices is observed for tq-Fds and -FD*/-Pbsj and also for the charm quark 
mass definitions tqMc^^'^ and r^Mc^^^ ; in the latter case, however, only when using ca by 
LANL. In general, on the coarsest two lattices we see significant discrepancies between the 
results of our non-perturbative computation obtained with the ALPHA-ca and the LANL- 
CA, respectively, which is in line with perturbative criterion on oMh consulted above. 

For the final results in the continuum limit of quenched QCD we quote the linear 
extrapolation of the data at the three finest lattice spacings, where we employ the ca 
of [20] for the improvement of the axial current. In this way we are well within the range of 
quark masses aMh, for which 0(a) improvement and thus the validity of the a'^-expansion 
and the associated scaling behaviour are expected to hold. In addition, the ambiguity 
owing to the choice of ca we observe has no significant effect at these lattice spacings. 

The fact that the continuum limit of the three definitions of the RGI charm quark mass 
(2.12) - (2.14) yield compatible continuum extrapolated results demonstrates nicely the 
universality of the continuum limit and thereby provides strong evidence for the correctness 
of the extrapolation procedure. Our final results in units of the hadronic scale tq are 
collected in the bottom row of table 4. The quoted errors include all errors, i.e. the 
statistical ones, those stemming from the fits at intermediate stages of the analysis as well 
as the uncertainties of ro/a and the various Z-factors. Using tq = 0.5 fm, we also convert 
our numbers to physical units: Fd, = 220(4)(?) MeV, mD|-mD, = 137(6)(?) MeV, M'i''"^ = 
1.59(2)(?) GeV, Mi^'^ = 1.60(2)(?) GeV and m]"^ = 1.61(2)(?) GeV. With the second error 
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we indicate the unknown systematic contribution due to the quenched approximation, by 
which also the result for the dimensionless ratio F^* / F^^ is affected. Note that thanks 
to the additional, finer lattice spacing and the substantially increased total statistics in 
the pseudo-scalar channel, we have gained a factor of about two in accuracy for vqMc and 
foPDs compared to the earlier investigations of [1, 14]. 

Whereas in the case of the RGI charm quark mass we observe agreement within errors 
with the result of ref. [14], the continuum limit of tq-Fds quoted in [1], which is based on data 
set A for the range of lattice spacings a « 0.08 — 0.05 fm, now proves to overestimate our 
result in table 4 for the leptonic decay constant of the Dg-meson by about three standard 
deviations. We have identified two reasons for this: Firstly, after increasing the statistics 
for ro-Fbs by generating data set B, the central value at a = 0.048 fm (/3 = 6.45) moved 
down by about 1.5o". Secondly, in addition to [1] wc now have results for a finer lattice 
spacing with a = 0.031 fm (/3 = 6.7859). Thus, increasing the statistics by combining data 
set A and B, discarding the data point at a = 0.079 fm (/3 = 6.2) from the continuum 
extrapolation and including the results at a = 0.031 fm instead causes a shift in the central 
value for F^^ by about three standard deviations compared to the result in [1]. Statistical 
fluctuations in conjunction with leading plus possibly higher-order lattice artefacts can 
misguide the true continuum extrapolation. 

This clearly reveals that it is indispensable to incorporate as small lattice spacings as 
possible into computations in the charm sector of 0(a) improved lattice QCD, in order to 
get a controlled handle on the continuum limit. 

5. Conclusion 

Full QCD simulations with improved Wilson fermions in large volumes and at small lattice 
resolutions very close to the physical point are now feasible [36,37]. Since their formulation 
is free of conceptual problems and since they are comparably cheap to simulate, they are 
now being applied on a large scale in the computation of Standard Model parameters at very 
high precision. The results of this work, although referring to the quenched approximation, 
emphasize the importance of performing a scaling study down to very small lattice spacings 
when calculating observablcs for hadrons containing a charm quark. Without the data 
at the finest lattice spacing and its satisfactory statistical accuracy obtained here, the 
continuum extrapolation would have yielded too large a value for Fd^ and also for Fd* / Fb^ . 

It is expected that these findings for the quenched theory carry over to the theory 
with dynamical sea quarks. We stress that only a scaling study of charmed observables 
down to extremely fine lattice spacings will allow for a reliable assessment of the cutoff 
effects. These findings, however, do not directly apply to other fermion discretizations and 
therefore, similar studies are crucial in each individual case. 
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